Characteristic Angular Scales in Cosmic Microwave Background Radiation 
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We investigate the stochasticity in temperature fluctuations in the cosmic microwave background 
(CMB) radiation data from Wilkinson Microwave Anisotropy Probe. We show that the angular fluc- 
tuations of the temperature is a Markov process with a Markov angular scale, Oiviarkov = LOIIq'^. 
We characterize the complexity of the CMB fluctuations by means of a Fokker-Planck or Langevin 
equation and measure the associated Kramers-Moyal coefficients for the fluctuating temperature 
field T(h) and its increment, AT = T(hi) — T(nz). Through this method we show that temperature 
fluctuations in the CMB has fat tails compared to a Gaussian distribution. 
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I. INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP) 
mission is one of the main experiments in the cosmic mi- 
crowave background (CMB) that determines the power 
spectrum of the temperature fluctuations in the CMB 
with high accuracy. It has been shown that the universe 
is geometrically flat, and that the dominant content of 
the universe is an exotic dark energy which causes the 
expansion of the universe to be accelerated [1-4]. The 
statistical properties of the CMB radiation is an im- 
portant tool for identifying the appropriate cosmological 
model and determining the parameters of the standard 
model [5] . The traditional way of studying the CMB data 
is through analyzing their angular power spectrum and 
computing the two-point correlation function. In order 
to gain full statistical information from the temperature 
fluctuations on the last scattering surface, we should use 
the n-point joint probability density function (PDF) or 
the n-point correlation function. 

At the same time, Gaussianity of the CMB anisotropy 
is an important question which is linked with the the- 
oretical predictions made by the inflationary cosmology 
[6-8]. With current developments in experimental CMB 
physics, we are now in a position to analyze very large 
data sets that provide information about large patches 
of the sky, measured with very high resolution and sensi- 
tivity. This means that we can precisely test whether 
the CMB is Gaussian. Tantalizing evidence of non- 
Gaussianity has been emerging in the WMAP sky maps, 
using a variety of methods. The statistical properties 
of the primordial fluctuations generated by the inflation- 
ary cosmology are closely related to the anisotropy of 



the CMB radiation. Thus, measurement of the possible 
non-Gaussianity of the CMB data is a direct test of the 
inflation paradigm. If the CMB anisotropy is Gaussian, 
then the angular power spectrum, or the two-point cor- 
relation function, fully specifies its statistical properties. 

In this paper we study the statistical properties of the 
CMB anisotropy through the n-point joint PDF of tem- 
perature fluctuations in the CMB data. We show that 
the fluctuations constitute a Markovian process with a 
Markov angular scale of OMarkov = l-Oljloo?) or m the 
Legendre space, ^Markov = at If confidence 

level. Using the Markov properties of the CMB data, a 
master equation for the angular evolution of the prob- 
ability distribution function - the Fokker-Planck (FP) 
equation - is obtained. To derive the relation between 
the standard power-spectral analysis and the Markov 
properties of the CMB data, we consider the joint prob- 
ability distribution P(T2,n2',Ti,ni) that describes the 
probability of finding simultaneously T\ in the direction, 
hi = (6*i , 0i ) , and T2 in the direction hz = (^2,^2)- We 
then evaluate the corresponding Kramers-Moyal (KM) 
coefficients, and compute the first and second KM coef- 
ficients (i.e., the drift and diffusion coefficients in the FP 
equation). We show that the higher-order coefficients in 
the KM expansion are very small and can be ignored. In 
addition, using the FP equation for the PDF, we show 
that temperature fluctuations on the last scattering sur- 
face are consistent with a Gaussian distribution. We also 
use the same analysis to investigate temperature fluctu- 
ations in the northern and southern hemispheres. The 
analysis detects some cold and hot spots in the southern 
and northern hemisphere, respectively. 

The organization of this paper is as follows: In Sec- 




FIG. 1. The arbitrary directions no, hi and hi in the sky 
to the last scattering surface. 

tion II we determine the 2-point joint PDF of tempera- 
ture fluctuations and obtain the correlation angular scale. 
In Section III we introduce the mathematical property 
of the PDF for the Markovian process and obtain the 
Markov angular scale for the WMAP data. We also de- 
velop the master equation in the form of a FP equa- 
tion for the CMB data. In Section IV we investigate the 
(non)-Gaussian nature of the temperature fluctuations in 
the CMB. The relation between the Markov and the an- 
gular power spectrum for this set of data is obtained in 
Section V. The conclusions are presented in Section VI. 



II. CORRELATION ANGULAR SCALE BY THE 
JOINT PDF DECOMPOSITION ANALYSIS 

The WMAP instrument is composed of 10 differencing 
assemblies, spanning five frequencies from 23 to 94 GHz 
[9]. The two lowest frequency bands (K and Ka) are 
primarily galactic foreground monitors, while the three 
highest (Q, V, and W) are primarily cosmological bands. 
The full width at half-maximum for the detectors is a 
function of the frequency and ranges from 0.82° at 23 
GHz to 0.21° at 94 GHz [10]. Here, we use the Inter- 
nal Linear Combination Map which is a weighted linear 
combination of the five WMAP frequency maps. The 
weights are computed using the criterion that minimize 
the galactic foreground contribution to the sky signal. 
The resultant map provides a low-contamination image 
of the CMB anisotropy. 

Let the temperature fluctuation in the CMB data in 
the direction of h = (0,<f) be represented by T(n), and 
define T(h) = (T(n) — T)/a, where T and a are the 
mean and variance of the temperature fluctuations, re- 
spectively. For a direction no with T(ho) = Tq, at any 



0.16 



0.15 - 



. 0.14 - 



CL, 
< 




0.13 - 



0.12 - 



FIG. 2. is the angular separation between two directions 
hi and 
hi. For Qc — 24.86° ± 0.03 we obtain the minimum value 
for the AP(T 2 ,Ti) = \P(T 2 , n 2 ; T u hi) - P(T U hi)P{Ti, h 2 )\ 

vs e. 



h ho, T(h) is given by a PDF P(T,h) (see figure 1). 
This means that T(n) for n ^ h is a stochastic (and 
possibly correlated) variable. 

Complete information about the stochastic process is 
obtained from the knowledge about all the possible k- 
point or, more precisely, fc-scale joint PDF (JPDF), 
P(Tk, hk] Tk-i, Hfe_i; Ti, hi), describing the probabil- 
ity of finding simultaneously, T\ in the direction hi, T2 
in the direction fi2, and so on up to 71- on the direction 
hk- Moreover, we can define the conditional probability 
distribution function. For a given T(ho) = Tq, the con- 
ditional probability of T{n) at the successive positions 
M i> f>-2, " ' ' , ft-hi to be T(hi), T{fi2), • • • , T{fik) is described 
by 



Pl(T k ,h k ; ■ ■ ■ ■ 1 T 1 ,hi\T ,ho)dT k dTk-i 

Prob{T(ni) e [T^T.+dT,] 

for i = l,2, ---k and T(h ) = T }. 



■dTi 



(1) 



In this notation, the superscript (for the P) represents 
the number of conditions (e.g., here, we have one condi- 
tion), while the subscript denotes the number of stochas- 
tic variables in the joint PDF. The fc-point JPDF of the 
temperature fluctuations is expressed by the product of 
the multiconditional PDFs as 

P (7fc_i,nfc-i;T fc _ 2 , rifc-2; Ti,hi;T , h ) = 
Pi ~ {Tk-x,fik-i\Tk-2, fik-2; T x ,hi;To, n ) 
x pfc- 2 (T , fe _2,n fe _ 2 |T fc _3,nfe-3; — ;T x ,ni;T ,ho) 



x...Pt(T 1 ,n 1 \T ,n )P(To,n ) 



(2) 



where, Pf(Ti, rai|To, no) denotes a conditional proba- 
bility of finding temperature fluctuation T\ in the di- 
rection hi, under the condition that the temperature 
T in the direction n has been found. This condi- 
tional probability can be expressed by 2-point JPDF, as: 
Pl(Ti,hi\T ,h ) = P( p(^fl°f o) , this result is impor- 
tant because, if the conditional PDF for the T\ is inde- 
pendent of To, then P\(Ti, ni|To, no) = P(T X , hi), in this 
case: P^Ti, n i; T , h ) = P(T U fn)P(T , h ). 

To derive the first characteristic angular scales in the 
CMB, we define the correlation angular scale, where the 
JPDF can be decomposed into the production of two in- 
dependent PDFs [11] as 



P{T 2 ,n 2 \Tx,hi)\e = P{Ti,h x )P{T 2 , 



»2, 



(3) 



Here, G is the angular separation between the two direc- 
tion hi and h 2 , and we assumed the statistical isotropy 
of the temperature fluctuations in the CMB in which the 
correlation function depends only on the angle between 
the two directions of hi and h2 [12]. In Fig. 2 we plot 
AP(T 2 ,Ti) = \P(T 2 ,n 2 ;Ti,hi) - P(Ti,hi)P(T 2 ,n 2 )\ in 
terms of @, using the least square method, 

2 f [P(T 2 ,h 2 ;Ti,hi) - PjT^h^PjT^h,)} 2 ^ 
X = 2 T — 2 a±id± 2 , 

J ® PDF + ^-joint 

where the minimum value of x 2 corresponds to the an- 
gular scale, at which the correlation disappears. Here, 
°~pdf an d o'i-joint are the variance of P(Ti,hi)P(T 2 , h 2 ) 
and P(T 2 , h 2 ;Ti,hi), respectively. Our analysis indicates 
that the minimum value of AP(T 2 ,Ti) happens at the 
e c = 24.86° ±0.03 with xl = 1-01 (xl = X /M, with TV 
being the number of degree of freedom). 



III. CMB DATA AS A MARKOV PROCESS: 
LEAST SQUARES TEST 

Now, let us check whether the CMB data follow a 
Markov chain and, if so, measure the second character- 
istic angular scale in the CMB, i.e., the Markov angular 
scale ©Markov - the scale over which the CMB data are 
Markov-correlated. In other words, the Markov angular 
scale OMarkov is the minimum angular separation over 
which the data can be represented by a Markov pro- 
cess [13-15]. The exact mathematical definition of the 
Markov process is given [16] by 

P(Tk,h k \T k _i,h k _i; ■ ■ ■ ; Ti, hi] T , n ) = 

P{T k ,h k \T k _i,n k _i). (4) 

Intuitively, the physical interpretation of a Markov pro- 
cess is that it "forgets its past," or, in other words, only 
the most nearby conditioning, say (T k , n k ), is relevant to 



the probability of finding a temperature T k at n k . In the 
Markov process the ability to predict the value of T(n) 
will not be enhanced by knowing its values in the steps 
prior to the most recent one. Therefore, an important 
simplification that can be made for a Markov process 
is that a conditional multivariate JPDF (Eq. 2) can be 
written in terms of the products of simple two-parameter 
conditional PDFs [16] as 



P(T k ,h k ;T k -i,n k - X ; ■ ■ ■ ;Ti,ni|Tb,n ) 

k 

~[P(T i ,h i \T i _i,h i _i). 

i=l 



(5) 



In what follows, we use the least-square method to 
determine the Markov angular scale of the CMB temper- 
ature fluctuations. Testing Eq. (4) for large values of k is 
out of our computational capability; however, for k = 3, 
where we have three vectors pointing the 2D celestial 
sphere, the Markov condition is as follows, 



P(T 3 ,h 3 \T 2 ,h 2 ;Ti,hi) = P(T 3 ,h 3 \T 2 ,h 2 ) 



(0) 



where Gy = arccos(ni.nj) is the separation angle of hi 
and hj and 631 > 032,021 (see Fig. 1). For simplicity, 
we let 02i = 032- A process is Markovian if Eq. (6) is 
satisfied for a certain angular separation which, in our 
notation, is ©32. We refer to it as the Markov angular 

SCale ©Markov)- 

In order to determine the Markov angular scale, we 
compare the three-point PDF with that obtained based 
on the Markov process. The three-point PDF, in terms 
of conditional probability functions, is given by 

P(T 3 , n 3 ; T 2 , n 2 ; Ti, hi) = 

P(T3, h 3 \T 2 , n 2 ; T 1; n 1 )T(T 2 , n 2 ; Ti, hi). (7) 

Using the properties of the Markov process and substi- 
tuting Eq. (6) we obtain: 



TMar(T3,n 3 ;T2,n2;Ti,ni) = 

P(T 3l h 3 \T 2 , h 2 )P(T 2 , n 2 ; Ti, h x ) 



(8) 



In order to check the condition for the data being a 
Markov process, we must compute the three-point JPDF 
through Eq. (7) and compare the result with Eq. (8). 
The first step in this direction is to determine the quality 
of the fit through the least-squared fitting quantity x 2 
defined by 

X 2 = J dT 3 dT 2 dTi [P(T 3 ,h 3 ;T 2 ,h 2 ;Ti,hi) 

~PM ar (T 3 ,h 3 ]T 2 ,h 2 ;Ti,hi)} 2 / [o-3-joint + ^Mar] ( 9 ) 
where CT I-joint 

and cr Mar are the variances of P(T 3 , h 3 ;T 2 ,h 2 ;Ti,hi) and 
-PMar(T3, h 3 ; T 2 , h 2 ; Ti, hi) , respectively. To compute the 
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FIG. 3. Xv in terms of angular scale O. The minimum 
value of xi — 1-388 corresponds to 0Markov = l.Oll ,; ,®, with 
la confidence level. 

Markov angular scale, we use the Likelihood statistical 
analysis [17]. In the absence of a prior constraint, the 
probability of the set of three-points JPDF is given by a 
product of Gaussian functions: 



p(arccos(rt 3 , hi)) = 



1 



FIG. 4. Relative likelihood function of the Markov angular 
scale of the CMB fluctuations, as a function of 0°. The in- 
tersections of the curve with the horizontal solid and dashed 
lines give the bounds with la and 2a confidence levels, re- 
spectively. The minimum value of Xv = 1-38 corresponds to 

©Markov = 1.011°;™, With la and ©Markov = 

2a confidence level. 
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[P(T3,n 3 ;T 2 ,h 2 ;Ti,ni) - PMm (T3,h 3 ;T 2 ,h 2 ;T 1 ,hj 



2 ( CT 3-joi„t + <ar) 



for three levels, 0.008 (red symbol), 0.005 (blue symbol), 
0.002 (green symbol). The conditional PDF P(T 3 |Ti), 
for T\ = ±<j are shown in the lower panel. All the scales 
are measured in unit of the standard deviation of the 
^temperature fluctuations. The CK equation, formulated 
— jn differential form, yields the following Kramers-Moyal 
(KM) expansion [16], 



(10) 



This probability distribution must be normalized. Evi- 
dently, when, for a set of values of the parameters, the \ 2 
is minimum the probability is maximum. Figure 3 shows 
the normalized xt as a function of the angular length 
scale = 032 = arccos(n3.n2). The minimum value of 
xl is 1.38, corresponding to 0Markov = l.Olj^o? w ^h Ic 
confidence level. Figure 4 shows the likelihood function 
of Markov angular scale of CMB. 

The Markov nature of the CMB enables us to derive 
a master equation - a FP equation - for the evolution of 
the PDF P(T, h), in terms of, for example, the direction 
n. One writes Eq. (8) as an integral equation, which is 
well-known as the Chapman-Kolmogorov (CK) equation: 

P(T 3 ,n 3 \Ti,hi) = J dT 2 P(T 3 ,h 3 \T 2 ,h 2 ) P(T 2 , n 2 |Ii, hi) 

(11) 

We checked the validity of the CK equation for de- 
scribing the angular separation of h\ and h 2 being equal 
to the Markov angular scale. This is shown in Fig. 5. 
In the upper panel we show the identification of the left 
(filled symbol) and right (open symbol) sides of Eq. (11) 



d 



^(-— )«[ J D(")(T,0)P(T,0)], (12) 



where i)W (T, (f>) are called as the KM coefficients. These 
coefficients can be estimated directly from the moments 
(M^ n >) and the conditional probability distributions as: 

DW(T,6) = ^7 lim M {n \ 
n! a</,^o 

M {n) = -^J dT '( T ' - T) n P(T', + A0|T, (/)). (13) 

According to the Pawula's theorem [16], for a process 
with L>( 4) - 0, all the D (,l) with n > 3 vanish, and the 
KM expansion reduces to the FP equation, known also 
as the Kolomogrov equation [16] 

d 2 



d 



P(T,, 



d_ 

df 



DW(T, 



Q T 2 



d<- 2 \t,, 



P(T,, 



(14) 

Here D^(T,(p) is the 'drift' term, describing the deter- 
ministic part of the process, while D^ 2 \T, cj>) is the 'diffu- 
sion' term. Using the Ito interpretation, the FP equation 
is equivalent to the following Langevin equation [16] 




FIG. 5. Verification of the validity of the Chapman-Kol- 
mogorov equation (Eq. 11) for the angular separation of ji\ 
and fi2 being equal to the Markov angular scale. In the upper 
panel we show the identification of left (filled symbol) and 
right (open symbol) sides of Eq. (11) for three levels, 0.008 
(red symbol), 0.005 (blue symbol), 0.002 (green symbol). The 
conditional PDFs P(T 3 |Ti), for Ti = ±cr are shown at the 
lower panel. All the scales are measured in unit of the stan- 
dard deviation of the temperature fluctuations. 

-^T(J>) = D^(T) + y/DW(T) /(</>). (15) 

Here, f(<fi) is a 5-correlated Gaussian random force with 
zero mean (i.e., < /(</>)/(</>') >= 3D5(</> - </>')). For the 
WMAP data, the drift and diffusion coefficients and 
£>( 2 ) are shown in Fig. 6. It turns out that the drift coef- 
ficient is a linear function in T, whereas the diffusion 
coefficient D^ 2 ' is a quadratic function. For large values 
of T, our estimates become poor and, thus, the uncer- 
tainty increases. From the analysis of the data set, we 
obtain the following approximate relations, 



£>U)(T) 



-0.330T, 



D {2) (T) = 0.060T 2 + 0.002T + 0.270, 



(16) 



where we used the isotropy assumption which implies 
that £>( n )(T,</>) = L> ( ")(T). The temperature field 
is measured in units of its standard deviation. The 
fourth-order coefficient D 1 - 4 ^ is, in our analysis, ~ 




FIG. 6. The drift and diffusion coefficients, D (1) (T) and 
D^ 2 \T) in Eq. (13), following, respectively, linear and 
quadratic behaviors. 

10~ 2 D( 2 \ so that we can ignore the coefficients D^ n > for 
n > 3. Furthermore, using Eq. (15), it becomes clear 
that we are able to separate the deterministic and the 
noisy components of the CMB fluctuations in terms of 
the coefficients and D^ 2 \ 



IV. THE (NON)-GAUSSIANITY TEST OF CMB 

Cosmological models of the structure formation are 
based on the standard inflationary paradigm, which im- 
plies Gaussian fluctuation of the density field [6-8,18]. 
The Gaussianity in the density perturbation directly 
translates into the Gaussianity of the CMB tempera- 
ture fluctuations. However, along with the standard in- 
flationary models, there exist theories that predict the 
non-Gaussianity of the primordial fluctuations. Inflation 
with two or more scalar fields can provide significant devi- 



ations from the Gaussianity [19-22]. Another possibility 
is the manipulation of the CMB data after the recom- 
bination, due to subsequent weak gravitational lensing 
[23,24] and various foregrounds, such as dust emission, 
synchrotron radiation, or unresolved point sources [25]. 
One should also take into account the additional instru- 
mental noise in the observational data [26]. 

There are standard methods of searching for non- 
Gaussian signature in the CMB data, such as using peak 
distributions [18,27], the genus curve [28,29], peak cor- 
relations [30] and global Minkowski functional methods 
[31]. More recent methods, such as a technique for 
studying hydrodynamic turbulence and detecting non- 
Gaussianity, and fractal analysis, can also be used for 
the CMB data [32,33]. Moreover, the Gaussianity of the 
CMB at different angular scales have been tested [34-57]. 
Most of the previous works were based on the consistency 
between the CMB data and simulated Gaussian realiza- 
tions. So far, they have found no significant evidence for 
cosmological non-Gaussianity. 

We note that for a Gaussian distribution, all the even 
moments are related to the second one through (T 2n ) = 
J^t(T 2 )" (e.g., for n = 2, (T 4 ) = 3(T 2 ) 2 ), while the 
odd moments are zero. We checked directly the relation 
between the moments for the CMB data. The results 
are summarized in Table I. The ratios (T 4 )/(3(T 2 ) 2 ) 
are, 1.077 ± 0.004 , 1.085 ± 0.005 and 1.048 ± 0.005, for 
the entire data set, and the northern and southern hemi- 
spheres, respectively. The odd moments (T 3 ) and (T 5 ) 
are not zero for different parts of the CMB data (see 
Table I). The odd moments of the temperature fluc- 
tuations can also be used as the indicators for the hot 
and cold spots. The southern part has a negative third 
moment, indicating greater prevalence of the cold spots, 
whereas the northern hemisphere has a positive third mo- 
ment, indicating domination of the hot spots. The ra- 
tio of the higher moments, such as (T 6 )/(15(T 2 ) 3 ), arc, 
1.398±0.026, 1.396±0.022 and 1.313±0.035 (for a Gaus- 
sian distribution this ration is unity) for the entire data 
set, and the northern and southern hemispheres, respec- 
tively. 

The ratio of the higher moments, 
such as (T 6 )/(15(T 2 ) 3 ), for the entire sky, and the north- 
ern and southern hemispheres are larger than those for a 
Gaussian distribution, implying that the probability dis- 
tribution function for the temperature fluctuations in the 
CMB has fat tails compared to a Gaussian distribution. 

Let us examine the predictions for the moments of the 
temperature fluctuations via the FP equation, and com- 
pare their values with the direct evaluation of results pre- 
sented in Table I. Using the stochastic properties of the 
temperature field, one can check the deviation of the fluc- 
tuations from the Gaussian distribution in the northern 
and southern hemispheres, as well as in the entire data. 
As pointed out in section III, the PDF of the temperature 
fluctuations satisfies the KM expansion. Using Eq. (12) 



we calculate the nth moment of the temperature fluctu- 
ations by multiplying both sides of Eq. (12) by T n and 
integrating over the temperature: 

ttCH = V i U ' „ {D (m) {T)T n - m ). (17) 

m—l v ' 

We set n = 4 in the above equation and find the equation 
for the fourth moment to be 

— IT 4 } = 4(L>«(T)T 3 ) + \2{D^\T)T 2 ) 
dtp 

+24(£> (3) (T)T) + 24(D (4) (T)). (18) 

The KM coefficients for the northern and the southern 
hemispheres are given in Table II. For the isotropic case, 
all the moments of temperature fluctuations are indepen- 
dent of 0, and the left-hand side of Eq. (18) vanishes. 
Using the results presented in Table II, for the northern 
and the southern hemispheres, Eq. (18) reduces to 

{T^north = (2.71 ± 0.06)(T 2 ) 2 + (0.028 ± 0.010)(T 3 ) y/J^), 
{T A ) SO uth = (3.27±0.09)(T 2 ) 2 - (0.017 ± 0.020) {T 3 )^ 2 ), 

(19) 

where for the northern and southern hemispheres, (T 3 ) = 
(1.283 ±0.125) x 10~ 5 , ^T 2 ) = (7.010 ±0.004) x 10~ 2 
and (T 3 ) = (-5.408 ±0.175) x lO" 5 , ^T 2 ) = (7.914 ± 
0.004) x 10~ 2 , respectively. We repeated the same pro- 
cedure for the entire data set and obtained the following 
relation between the moments 

(T 4 ) = (3.07 ± 0.07) (T 2 ) 2 + (0.034 ± 0.020) (T 3 ) ^AT 1 ), (20) 

where (T 3 ) = (-2.032 ± 0.108) x 10~ 5 and ^/JT 2 } = 
(7.475±0.003) x 10" 2 . The drift and diffusion coefficients 
of the northern and the southern hemisphere arc shown 
in Fig. 7. 

As mentioned earlier, for a Gaussian distribution, all 
the even moments are related to the second one, while 
the odd moments are zero. Eqs. (19) and (20) show that 
we have deviations from the Gaussianity, resulting from 
the coefficient of (T 2 ) 2 . Note that the result presented 
by Eq. (20) are compatible with direct calculation of 
moments for different part of the CMB data. 

V. MARKOV NATURE OF CMB AND 
ANGULAR POWER SPECTRUM 

Using statistical isotropy, the two-point temperature 
correlation function is expanded in term of the Legendre 
functions as 

(T( ni )T(n 2 )) = -±-C,Pi(coB(e)), (21) 
i 
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FIG. 7. Drift and diffusion coefficients D (1) (T) and 
£)( 2 ) (J 1 ) f or tj le northern and the southern hemisphere parts 
of the CMB data. 



where C/ is the angular power spectrum, Pi is the Leg- 
endre polynomial of order I and 8 = arccos(ni.n2)- To 
find the relation between the Markovian properties of the 
CMB data and the Cj's, we start with the temperature- 
increment moments, 52 = ([T(ni) — T(fi2)] 2 ). The sec- 
ond moment S2 allows us to determine (T (ni)T (712)) and, 
hence, the Cj's. On the other hand, 5*2 can be obtained 
through P(AT, 8), where AT = T(n x ) - T{h 2 ). The FP 
equation [16] for AT(6) is given by, 
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FIG. 8. Comparison of the numerical solution of Eq. (24) 
and direct computation of (T (rii)T (712)) ■ 



T> (2) (AT,8) = [0-021 + 0.025 exp(— 



e 



5.896 



)](AT) 5 



+0.279 



0.014 

Q0.429 ' 



(23) 



Here, AT is measured in units of the standard deviation 
of T. Equation (22) allows us to obtain an equation for 
the second moment of the temperature increments, S2 ■ 
Multiplying Eq. (22) by AT 2 and integrating over the 
increment of the temperature yield 



— <T(ni)T(n 2 )} 



2[a(e)+w(e)KT(ni)T(n a )) 



[2a(8) + 2^(8) -/3(8)], 



(24) 



where a 2 = (T 2 ), DW(AT,0) = a(8)AT(8), and 
(AT, 8) = /3(8)cr 2 + A(8)aAT(8) + w(0) AT(8) 2 . 
The coefficients a(Q), /3(8), A(8), and w(8) are given 
by Eq. (23). 

We obtain the correlation functions through numerical 
solution of Eq. (24) and compare it with the direct corre- 
lation shown in Fig. 8. The two approaches yield similar 
behavior. At 8 ~ 39° the correlation function vanishes. 
The relation between the KM coefficients with Ci is then 
obtained by expanding the solution of Eq. (24) in terms 
of P;. 



P(AT, 8) 







d AT 



T>W(AT,8) 



if 



9 AT 2 



D (2) (AT,e)]P(AT,8). (22) 



From the analysis of the CMB data we obtain the follow- 
ing equations for £>W(AT, 8) and D^(AT, 8) 



D (1) (AT,8) = (-0.190- 



0.182 
8 



) AT, 



VI. CONCLUSION 

We studied the stochastic nature of the temperature 
fluctuations in the CMB. The Markov angle, as the char- 
acteristic scale of the Markov properties of the CMB 
data, was obtained. According to the theory of stochas- 
tic process, the CMB data at scales larger than the 



Markov angle can be considered as a Markov process. 
This means that the data located at separations larger 
than the Markov scale can be described as a Markov 
chain. We also obtained the angular scale correlation 
for the CMB data and showed that it is independent of 
the Markov angular scale. This point can be explained as 
a simple example of the Brownian motion. Its dynamics 
is described by the Langevin equation, 4| = —av + 
(similar to Eq. (15)), where the force rj(t) is a zero-mean 
Gaussian white noise. It is known that this process is a 
Markov process with a Markov time scale of unity [16]. 
On the other hand, the correlation time scale is of the 
order of 1/a, which is the characteristic time scale in the 
Langevin equation [16]. This means that for a Marko- 
vian process the Markov time scale does not depend on 
the correlation time scale. 

We showed that the probability density of the tem- 
perature increments satisfies a Fokker-Planck equation 
that encodes the Markov property of the fluctuations. 
We gave the expressions for the Kramers-Moyal coeffi- 
cients of the stochastic process T and AT(9), using the 
polynomial anzats [58-63]. The Fokker-Planck equation 
enables us to derive a simple equation that governs the 
phenomenon in terms of azimuthal spherical coordinate. 
One of the important points that can be tested by this 
method is the (non)-Gaussianity of the temperature field. 

Using the Fokker-Planck equation, we obtained the re- 
lation between the fourth and the lower moments, and 
observed small deviations from the Gaussianity. The 
same calculation was carried out by dividing the data 
into the northern and southern hemispheres. The PDF 
of the temperature exhibits fat tails for the northern 
and the southern hemispheres. The third moment in- 
dicates that, we have hot spots in the north, in contrast 
with the cold spots in the south hemisphere. To show 
a link between our method and the standard analysis of 
the CMB data through C; calculation, we obtained the 
evolution equation for the correlation function through 
the Fokker-Planck equation that governs the tempera- 
ture increments. We found good agreement between our 
method and that of direct correlation function calcula- 
tion. 

The most important result of this paper is a possible 
interpretation of the CMB data in terms of the Markov 
angular scale. In this work we computed the Markov an- 
gular scale for the temperature fluctuations in the CMB, 
in the range of LOT^oo?' where comparison with the 
event horizon shows that the two angles are of the same 
order of magnitude. A possible interpretation may be 
the physical connection of the Markov angular scale for 
the CMB data to the event horizon at the last scattering 
surface. 

From the definition of the event horizon, two points 
located further apart than the Markov angular scale do 
not connect gravitationally. From the inflationary sce- 
nario for the early universe, all the points in the universe 



have been correlated during the inflationary epoch, and 
get almost uniform Harrison-Zel'dovich spectrum. Af- 
ter the end of the inflation, the perturbation remained 
from the inflation started to grow and the gravitational 
effect of each point could travel within the event horizon 
scales, so that outside the horizon we have the primor- 
dial spectrum from the inflation, while inside the horizon 
this spectrum has been vanished due to the gravitational 
interaction between the points. The result is that, if we 
look at the density perturbation in the CMB as a stochas- 
tic field, the conditional probability of finding two points 
inside the event horizon scale will depend on all the par- 
ticles in between, while for the points outside the event 
horizon, the memory of scale-invariant spectrum should 
be restored. A better interpretation of this issue needs 
a A r -body simulation of structure formation from the in- 
flationary epoch to the last scattering surface, together 
with a simultaneous calculation of the Markov length. 
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TABLE I. The values of Moments (T n ) for whole sphere, northern and the southern hemisphere of the CMB map data. 





(T 2 ) 


(T 3 } 


(T 4 ) 


(t b } 


Whole-sphere 


(5.588 ± 0.005) x 10^ 


(-2.033 ±0.108) x 10~ s 


(1.009 ±0.003) x 10~ 4 


(-1.298 ±0.132) x 10"° 


Northern- hemisphere 


(4.914 ± 0.006) x 10" 3 


(1.283 ±0.125) x 10" b 


(7.867 ± 0.0305) x 10^ 


(7.632 ± 1.028) x 10"' 


Southern-hemisphere 


(6.264 ± 0.007) x 10" 3 


(-5.408 ±0.175) x 10 _& 


(1.233 ± 0.005) x 10" 4 


(-3.392 ± 0.244) x 10" b 



TABLE II. The values of Kramers-Moyal coefficients for the northern and the southern hemisphere of the CMB map data. 







L» (2) (T) 


£)(3)( T ) 


D (4) (T) 


Northern- hemisphere 


-0.370T 


0.290-0.003T±0.075T 2 


-0.001-0. HOT 
±0.003T 2 -0.010T 3 


0.040-0.007T±0.0182T 2 
-0.0005T 3 ±0.001T 4 


Southern-hemisphere 


-0.290T 


0.250±0.007T±0.047T 2 


0.004-0.078T 
-0.005T 2 -0.003T 3 


0.031-0.003T±0.020T 2 
±0.001T 3 -0.001T 4 



